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Abstract 

We study the electron transport through a graphene nanoribbon-superconductor junction. Both 
zigzag and armchair edge graphene nanoribbons are considered, and the effects of the magnetic 
field and disorder on the transport property are investigated. By using the tight-binding model and 
the non-equilibrium Green's function method, the expressions of the current, conductance, normal 
tunneling coefficient, and Andreev reflection coefficient are obtained. For a clean system and at 
zero magnetic field, the linear conductance increases approximatively in a linear fashion with the 
on-site energy. In the presence of a magnetic field and a moderate disorder, the linear conductance 
exhibits plateau structures for both armchair and zigzag edges. The plateau values increase with 
the width of the graphene ribbon. With a wide sample width, a saturated plateau value of \v\e^ /h 
emerges at the filling factor u. For a small filling factor, the conductance can reach the saturated 
value at a small width, but for a high filling factor, it requires to have a quite wide sample width 
to reach the saturated value. In particular, the Andreev reflection coefficient is always at 0.5 after 
reaching the saturated value, independent of any system parameters. In addition, we also consider 
the finite bias case, in which the Andreev reflection coefficient and normal tunneling coefficient are 
studied. 



I. INTRODUCTION 



The recent experimental realization of graphene,-*^ a single-layer carbon atoms ar- 
ranged in a honeycomb lattice, has generated a great attention in the condensed matter 
community.-!^ Graphene has an unique band structure with a linear dispersion relation 
of the low-lying excitations, which leads to many peculiar properties,-!^ such as its quasi- 
particles obeying the Dirac-like equation and having the relativistic-like behaviors with zero 
rest mass, the Hall plateaus having the half-integer values g{n + l/2)e^ /h with the degen- 
eracy (7 = 4. For the neutral graphene, its Fermi level passes through the Dirac points, the 
six corners of the hexagonal first Brillouin zone. By varying the gate voltage, the charge 
carriers of graphene can be easily tuned experimentally. Then the Fermi level can be above 
or below the Dirac points. 

Very recently, some works have begun to investigate the transport behaviors of the 
graphene-superconductor junctions.-'^'^'^i^iii^i'iiii^'i^ii^'i^'i^'i^ While a metal coupled to a 
superconductor, the Andreev reflection occurs in the interface between the metal and 
superconductor,^^ in which the interface reflects an electron incident from the normal metal 
side as a hole and a Cooper pair is created in the superconductor. For a bias below 
the superconductor gap, the Andreev reflection determines the conductance of the metal- 
superconductor junction since the normal tunneling cannot occur. In the usual metal- 
superconductor junction, the Andreev reflected hole retraces the path of the incident elec- 
tron, so this Andreev reflection is also called Andreev retroreflection. But for the graphene- 
superconductor junction, Beenakker recently found that a new kind of reflection (specular 
Andreev reflection) occurs while the incident electron and reflected hole are at the conduction 
and valence bands, respectively.- Afterwards, many papers have studied the graphene and 
superconductor hybrid system, including the graphene-based normal-superconductor (N- 
S),i^;^iiiii S-N-S,-i^'i2 S-insulator-S,~ S-ferromagnet-S,ii and etc. Several other effects due 
to the coupling of graphene and superconductor, such as Josephson effect^"^*^ and multiple 
Andreev reflection processes,— have been theoretically analyzed. On the experimental side, 
good contacts between the superconductor electrodes and graphene have been realized,— 
and the Josephson current through a S-graphene-S junctions has been measured. A finite 
supercurrent was observed at zero charge density. 

In this paper, we carry out theoretical study of the transport characteristics of a graphene 
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nanoribbon-superconductor junction. There are four new aspects beyond the previous stud- 
ies: (i) We study the system consisting of a graphene nanoribbon with a finite width coupled 
to the superconductor electrode. The previous theoretical papers only consider the infinite- 
wide graphene-superconductor junction or a graphene strip between two superconductor 
leads with the strip width much larger than the strip length. On the experiment side, the 
graphene nanoribbon has been successfully fabricated, and the width of the nanoribbon can 
be in the order of ten or sub-ten nanometers;^^ (ii) In our model, the graphene nanorib- 
bon is directly coupled to the superconductor electrode, and the incident electrons from the 
graphene are allowed to enter into the superconductor electrode as the Cooper pairs. In the 
previous papers, those authors only considered a pair potential in the graphene induced by 
depositing of a superconductor electrode on top of the graphene sheet, (iii) We consider a 
perpendicular magnetic field applied to the graphene, as was done in a recent experiment.— 
On the superconductor side, the magnetic field vanishes due to the Meissner effect, (iv) The 
effect of disorder on the transport property is investigated since in a real graphene sample, 
the disorder is always there to a certain degree. In fact, in the previous studies, the effects 
of disorder and magnetic field are thus far neglected. 

By using the tight binding model and the nonequilibrium Green function method, the 
current and Andreev reflection coefficient are obtained. Both zigzag edge and armchair edge 
graphene nanoribbons are considered. For the zigzag edge and at a zero magnetic field, 
the linear conductance exhibits step structures for the narrow graphene ribbon. With a 
magnetic field, the conductance depends strongly on the system parameters. In the presence 
of disorder, the linear conductance shows plateaus at a high magnetic field. On the other 
hand, for the armchair edge, a zero-conductance region emerges because of the existence of 
an energy gap in the graphene nanoribbon. This zero conductance is robust against disorder. 
In addition, we also consider a finite-bias case, in which the Andreev reflection coefficient 
and normal tunneling coefficient are investigated. 

The rest of the paper is organized as follows. In Section II, the model for graphene 
nanoribbon-superconductor junction is presented and the formalisms for calculating the 
current and the Andreev reflection coefficient are derived. In Section III and IV, we study 
the linear conductance and the transport with a finite bias, respectively. Finally, a brief 
summary is presented in Section V. 
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II. MODEL AND FORMALISM 



We consider the system consisting of a graphene nanoribbon coupled to a superconductor 
lead (as shown in Fig.l) with the Hamiltonian: 

H = Hg + Hs + He, (1) 

where Hq, Hs, and He are the Hamiltonians of the graphene region, superconductor lead, 
and coupling of the graphene and superconductor lead, respectively. For a semi-infinite 
graphene nanoribbon, Hq in the tight-binding representation is of the form:— 

= I] eifllaia - te'^^'alaja, (2) 

i,u <ij>,(7 

where Ojo- and a]^ are the annihilation and creation operators at the discrete site i, and is 
the on-site energy which can be controlled by the gate voltage in an experiment. Two kinds 
of edges, zigzag and armchair, are considered (see Fig. la and lb). The graphene ribbon is 
divided into two regions. The left-side of the semi-infinite region is without disorder and 
ej = El there. The disorder exists only in the center region of the graphene-nanoribbon 
(see the box with the dotted-line in Fig.l). Here we consider the on-site disorder causing by 
the nonmagnetic impurities or by the random potential difference of the substrate. Due to 
the disorder, the on-site energy ei = + Wi, where Wi is the on-site disorder energy and 
Wi is uniformly distributed in the range [—W/2, W/2] with W being the disorder strength. 
The size of the disorder region is described by the width and length L. In Fig. la and lb, 

= 3, L = 4 and N = 4, L = 2, respectively. The second term in Eq.(2) is the nearest 
neighbor hopping. When the graphene ribbon is under a uniform perpendicular magnetic 
field B, a phase (pij is added in the hopping elements, and (pij = Ji A - dl/cpo with the vector 
potential A = {—By, 0,0) and 0o = ^/e. 

Experimentally, it is possible to have the superconductor electrode in a good contact with 
the graphene.— The electrons from the graphene can easily enter into the superconductor 
electrode as the Cooper pairs or vice versa. So we consider that the graphene nanoribbon is 
directly coupled to the superconductor electrode. The superconductor electrode is described 
by a continuum model and it does not have the honeycomb structure of the graphene. Then 
the Hamiltonian Hs is: 

Hs = Y: ekC&k. + E (A^kT^-ki + A&-ki&kT) , (3) 

k,cr k 
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where feko- and b]^^ are the annihilation and creation operators in the superconductor lead 
with the momentum k = {kx,ky). Here we consider the s-wave superconductor and A 
is the superconductor gap. The superconductor region is without the magnetic field due 
to the Meissner effect or that the magnetic field is only added in the graphene region. 
The Hamiltonian He of the coupling between the superconductor lead and the graphene 
nanoribbon is: 

Hc = Y.tca\X{yi) + H.c. (4) 

i,C! 

Here only the surface carbon atoms couple to the superconductor lead, and Ui is the vertical 
position of the carbon atom i. ba{y) is the annihilation operators at the position (0,?/) of 
real space, and 

ha{y) = E e^'^'^k. (5) 

The current flowing through the graphene nanoribbon-superconductor junction can be 
calculated from the evolution of the total number operator for electrons in the left graphene- 
nanoribbon lead,— 

d 



duo 

/ 

where tij = td^^^ . Here i E L and j E C represent that the site index i and j are in 
the left graphene lead and center region. Gfj{uj) is the matrix Green function in Nambu 
representation, and it is the Fourier transformation of G^(t): 



^(a]^(0)a,t(t)) (%i(0)a,T(t)) 



{a]^{0)al{t)) (a,i(0)aj^(t)) ) ' 
By using the Dyson equation, the current expression in Eq.(6) can be rewritten as: 



(7) 



/ = ^|d^Tr| ^ ](g)I^a(S2-I]l)G< + S<(G'--G'^)}. (8) 



Here G^'"'^(co') are the 2Nc x 2Nc matrix Green's functions in the center region with Nc 
being the number of sites in the center region. The retarded and advanced Green's functions 



Qr,a defined in the standard way.— Inc is the Nc x Nc unit matrix. are the 

retarded, advanced, and lesser self-energies of couphng to the left graphene lead, and they 
are: 

\ 




(9) 



^inSnm,ll^'rnj 

''inS>nm,22''mj J 

S2 = (10) 



(11) 



tf^{u)TLi{uj) ) 

where fi{uj) = f{uj—eV) and fi{uj) = f{uj+eV) with V being the bias voltage and /(cu) being 
the Fermi distribution function, and T^{u) = ^(S^ — S2)ii and Tli^^u) = i(S'£ — S2)22- 
Snm('^) Eq.(9) is the surface Green's function of the semi-infinite graphene nanoribbon, 
that can be numerically calculated.^ With the aid of the self-energy functions in Eqs.(9-ll), 
the current / is finally reduced to: 

I = Il + h. (12) 

if f 

/t = jj duTrTL^{G< + /t(G'- - G")}n, (13) 

h = - jj dujTrTLi{G< + /^(G'^ - G'^)^. (14) 

As shown in the Appendix, the self-energies S^"^'*^ of coupling to the superconductor lead 
can be obtained by = —{i/2)TR, = (i/2)r^, and = if{uj)Tji. Then by using 
the Keldysh equation G< = G"S<G", G"^ - G'^ = G''(S" - S'^)G^ and the self-energies 
5]r,a,< _ 5]2"'^ + S^"'^, the currents /| and can be rewritten as: 

h = 11 duTr{T,^ [G^TaG\, (f) - f) + Tl^GI^T^^GM ' /i)}, (15) 
h = -{j du:Tr{TLi [G^T^G^^^ Ui - /) + r^iG^ir^TGi2(/i - /t)}- (16) 

Here Tr{TL^ [G^TrG\-^} = T^{uj) and Tr{TLi [G'^TrGX^} = ^i(^) are the normal tun- 
neling coefficients for the incident spin-up electron and spin-down hole with the energy 
u, and TriTL^Gl^TLiG^i} = Ta^{uj) and Tr{rLiGl^T^G^2} = Tai{uj) are the Andreev 
reflection coefficients. Since the Pauli matrices (Tx,j/,^ commute with the Hamiltonian H, 
the normal transmission coefficients T-^iuj) = T|(— cu) = T{uj) and the Andreev reflection 
coefficients Ta-\{uj) = T4|(— cu) = Ta{uj). 
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In the following, we need to calculate the Green's functions C and of the center 
region. Since the self-energy S'" has been obtained before and by using the Dyson's equation, 
the Green's function is simply of the form 

G'^(^) = 1/ {uhNc - Hcenter ' ^1 , (17) 

and in addition G" = G*"^, where 'Reenter is the Hamiltonian of the center region in the 
Nambu representation. 

In the numerical calculations, we take the hopping energy t = tc = 2.75eV and the 
nearest-neighbor carbon-carbon distance a = 0.142nm as in a real graphene sample.-i^ The 
superconductor gap A is set to A = t/2750 = ImeV, and the Fermi wave-vector kp = lA ^ . 
The temperature T is set to zero since T can be as low as IK in a real experiment and thus 
ksT is much smaller than all other relevant energies, such as t and A. The magnetic field is 
expressed in terms of (p with (p = (3v^/4)a^i?/0o and (3v^/2)a^-B is the magnetic flux in 
the honeycomb lattice. In the presence of disorder, the curves are averaged over up to 1000 
random configurations. 

III. THE LINEAR CONDUCTANCE 

In this section, we consider the small bias limit and investigate the linear conductance. 
When the bias V is smaller than the gap A, the normal tunneling processes can not occur 
and T{uj) = for Ic^l < A. Then only Andreev reflection processes contribute to the current, 
and the linear conductance G = hm^dl/dV = {Ae^ / h)TA{0) at zero temperature. In the 
following, we carry out numerical studies of graphene nanoribbons with both zigzag and 
armchair edges. 

A. The zigzag-edge case 

First, we study the clean graphene nanoribbons with the disorder strength W = 0. Fig. 2 
shows the linear conductance G versus the on-site energy (i.e. the energy at the Dirac 
point) with and without the magnetic field. The energy can be controlled by the gate 
voltage in an experiment. For El > 0, the charge carrier of graphene is hole-like, and it 
is electron-like for El < 0. In the absence of a magnetic field (0 = 0), the conductance 
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G is approximately linear with \El\ due to the linear increasing of the carrier density. For 
a narrow graphene nanoribbon (e.g. = 40), the conductance G clearly shows the step 
structures because of the sub-bands from the finite width. When a sub-band passes through 
the Fermi energy Ep {Bp = 0), a step appears. For = 40, the width of the graphene 
ribbon is (3iV — l)a 17nm. The graphene ribbon with this width has been realized in a 
recent experiment.— On the other hand, for a wide graphene nanoribbon (e.g. = 70), the 
step structures faint away due to the reduction of the interval of the sub-bands. 

While in the presence of a strong magnetic field, the conductance G does not show a clear 
pattern and depends strongly on El and the width (see Fig. 2b). Raising the disorder from 
zero, the conductance G in the small-value region is increased while G in the large-value 
region is decreased (as shown Fig.3), meanwhile some plateaus emerge at moderate disorder 
strength, e.g W = 2M These plateaus origin from the mixture of the electron and hole edge 
states, which will be discussed in detail in the last paragraph in this sub-section. 

Fig. 4 shows the linear conductance G versus the on-site energy Ep with a moderate 
disorder strength W. G exhibits the plateaus with or without a magnetic field. In absence 
of the magnetic field, the conductance G is similar to the disorder- free case (compare Fig. 2a 
and Fig. 4a), and the plateaus of the conductance are equal-spaced in energy. These plateaus 
are from the discrete sub-bands. In a graphene sample, due to the linear dispersion relation, 
the sub-bands are equal-spaced, so are the plateaus. For a wider graphene ribbon, the sub- 
bands are closer, then the width of the plateaus are smaller, so that the plateaus are fainted 
at large width (e.g. A^ = 70). On the other hand, in the presence of a magnetic field, the 
width of conductance plateaus are independent to the width A^ of the graphene ribbon, and 
the plateaus are always clear regardless of A^. Now the plateaus are equal-spaced in the 
scale of Ep, and the values of the conductance plateaus are determined by the filling factors 
u of the Landau level and the width A^ of the graphene ribbon. The wider A^ is, the larger 
the conductance value is. But the conductance reaches a saturated value \v\e^ /h at large A^ 
(see Fig. 5b). Fig.5 shows the conductance G versus the width A^ of the graphene ribbon. 
For = 0, G increases approximatively in a linear way with A^ (see Fig.Sa). But at high 
magnetic field G has a saturated value (see Fig.Sb). For small filling factor(e.g. Ep = O.lt 
with u = 2), G reaches the saturated value with small A^ (A^ = 40). For large filling factor, 
G reaches the saturated value only with quite large A^. 

With the aid of the edge states, these phenomena can be well explained. With a high 
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magnetic field, the edge states that carry charges are formed. In the interface of the graphene 
and the superconductor, the edge states extend from one boundary to the other along the 
interface, in which the Andreev refiection occurs. So the wider the graphene ribbon is, the 
larger the probability is for the Andreev reflection. In the large N limit, the electron and hole 
edge states are well mixed, thus, the Andreev reflection coefficient is 0.5, independent of any 
system parameters, such as the width N, the on-site energy El, the magnetic field strength 
0, and the disorder strength W. Then the conductance G — (2e^//i)2TA(0) = |//|e^//i. 

B. The armchair-edge case 

In this sub-section, the linear conductance G in the armchair edge case is numerically 
investigated. Fig.6 shows the conductance G versus the on-site energy El. Without a 
magnetic field (0 = 0), G increases linearly with the energy \El\ in the absence of disorder 
(see Fig. 6a). The disorder evidently enhances the conductance G in the small \El\ region 
(see Fig. 6b), thus, G departs from the hnear relation with \El\. In contrast with the zigzag 
edge case, it has two obvious characteristics: (i) There is a zero conductance G region near 
El = ioT N = 3m or 3m -|- 1, because that an energy gap emerges at the armchair edge 
graphene ribbon causing the Andreev reflection to vanish. This zero conductance G region 
still exists in the presence of disorder (e.g. W = 2). (ii) The step structures from the 
sub-bands are not apparent, although the width of the graphene ribbon is y/3Na ~ 17nm 
for N — 70. For the zigzag edge case with this width the step structures are clearly seen 
(see Fig. 2a and Fig.4a). 

With a magnetic field, the Landau levels are formed and the conductance G departs com- 
pletely from the linear relation with \El\. For the clean system {W = 0), the conductance 
is quite small at the smallest filling factor |z/| = 2, and exhibits some peaks at the higher 
filling factors — 6, 10, 14, etc (see Fig.6c). On the other hand, in the presence of disorder 
(W — 2), the conductance G shows plateaus and the plateau values are \v\e^ /h in the large 
width N limit. This is because of the mixture of the electron and hole edge states and the 
Andreev refiection coefficient is 0.5 at large A^. The characteristics of the plateaus at the 
moderate disorder strength are similar to that of the zigzag edge graphene ribbon. 

Fig. 7 shows the conductance G versus the width N of the graphene ribbon for a moderate 
disorder strength. At zero magnetic field, the conductance G increases linearly with the 
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width N as appears in a classical system. But at a high magnetic field, although the 
conductance G still increases with the width, a saturation value \v\e^ /h appears, same as in 
the zigzag edge case. 

IV. THE FINITE BIAS CASE 

In this section, the case with a finite bias is investigated. With the bias F > A, the 
normal tunneling processes also occur, and the current / is: 

7=^1 duj{nuj){f^ - /) + Ta{uj)U^ - /i)}. (18) 

Following, we numerically study the normal tunneling coefficient T{uj) and Andreev reflec- 
tion coefficient Ta{uj) for the zigzag edge graphene ribbon. Fig. 8 shows T{ijj) and Ta{oj) 
versus the energy uj of the incident electron for the clean system. The normal tunneling 
coefficient T{uj) is zero when \uj\ < A because of the superconductor gap, and T(co') is near 
1 at \uj\ > A since there is no barrier at the interface of the superconductor and graphene. 
Next, we focus on the Andreev refiection coefficient Ta{uj). For zero magnetic field with 
= 0, Ta{uj) is almost zero for \uj\ > \El\ (see Fig. 8a), implying that the specular An- 
dreev refiection is very weak at = 0. But the usual Andreev retroreflection still occurs, 
and Ta{uj) is quite large for < \El\. With a magnetic field (see Fig. 8c), both kinds 
of Andreev refiections occur simultaneously, and Ta{uj) is always finite regardless whether 
\uj\ < \El\ or \uj\ > \El\. Ta{uj) has a peak at = ±A and quickly decays for \uj\ > A, 
which is similar to a normal metal-superconductor junction.— 

Finally, the effect of the disorder on the normal tunneling coefficient T{uj) and Andreev 
refiection coefficient Ta{uj) are studied. The normal tunneling coefficient T{uj) is almost 
unaffected by a moderate disorder strength W, T{uj) is still zero for |ci;| < A and near 1 for 

> A (see Fig. 9b and 9d). However, the Andreev reffection coefficient Ta{uj) is evidently 
affected by the disorder (see Fig. 9a and 9c). Both specular Andreev reffection and the usual 
Andreev retroreffection occur, and Ta{uj) is close to 0.5 in the whole range of \uj\ < A. 

V. CONCLUSION 

In summary, by using the non-equilibrium Green's function method, the electron trans- 
port through the graphene nanoribbon-superconductor junction is investigated. Both zigzag 
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and armchair edge graphene nanoribbons are considered. The effects of a magnetic field and 
disorder on the transport property are discussed. In the clean system and without a mag- 
netic field, the linear conductance increases approximatively in a linear fashion with the 
on-site energy for the case with the armchair edge or the wide zigzag edge. In the pres- 
ence of a magnetic field and moderate disorder, the linear conductance exhibits the plateau 
structures for both armchair and zigzag edge nanoribbons. The plateau value increases with 
the width of the graphene ribbon, but reaches a saturation at |i/|e^//i (i^ is the filling factor) 
for the wide graphene ribbon. In addition, the case with a finite bias is studied, and the 
dependence of the Andreev reflection and normal tunneling coefficients on the energy of the 
incident electron are discussed. 
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APPENDIX 



In this appendix, wc derive the surface green function of the superconductor lead and 
the self-energies S^"'^ of coupling to the superconductor lead. The definition of the surface 
green function is: 

[{{bl{y,t),b\{y',0)}) {{b\{y,t)My'm 

and gg{y,y',u) is the Fourier transformation of gg{y,y',t), where y and y' are the real- 
space positions on the surface of half-infinite superconductor lead. Applying the equation 
of motion, g,s{y,y','^) can be written as: 



A a;+ 



11 



where ojj^ = u + iO+. Next, we calculate the sum, 
6k)/K - 4 - A') or F(k) = A/K - 4 - A'), 



k 



iky{y~y') 



dd / dkkpke'^^y-y"^'''^^F{k) 

IT J 

deuUk{y-y'))p{eu)F{k). 



where Jq is the Bessel function of the first kind, pk is the density of state in the k space 
and p(efc) = 2T^kpk{dk/ dek) is the density of state in the energy space. In the above steps, 
we have assumed the s-wave superconductor so that ek only depends on k = |k|. In the 
following, we assume that the density of state p(efc) = p is independent of the energy and 
Jo{k{y — y')) only depends on the Fermi wave- vector kp.— These assumptions are reasonable 
because that the main contribution to the transport behaviors is these electrons with their 
energies near the Fermi energy. Then gg{y,y' ,00) reduces to: 

/ 



Ssiy^y'^^) = Mkpiy - y'))p J dek 



1 



UJt, 



A2 



uj+ + ek A 
A uj> 



By using the theorem of residue, the integration / dek in the above equation can be obtained 
analytically,— and the surface Green's function gg{y,y',u) changes into: 

gs(i/> y ,^) = -i-^pJoikpiy - y 

\A/uj 1 

where (3{uj) = \uj\/^/lu'^ — A"^ while Ic^l > A and l3{uj) = lu / A"^ — u;^) while Ic^l < A. 
After solving the surface Green's function gg(y,y',uj), the self-energies 5]^'^'^ are obtained 
straightforwardly. 



= -mp\tc\^jQ{kF{yi-yj))l3{uj) 

= -(V2)r^,,,(cu), 



A/u 



\A/uj 1 



S^(cu) = it/2)Tn{uj) and S<(a;) = tfico)Tnico). 
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FIG. 1: (a) and (b) arc the schematic diagrams for the zigzag and armchair edge graphene 
nanoribbon-superconductor junctions, respectively. 

FIG. 2: (Color online) The linear conductance G vs. the energy for different width N at the 
clean system with W = Q. The panels (a) and (b) are for the magnetic field strength ^ = and 
= 0.007, respectively. 

FIG. 3: (Color online) The conductance G vs. El for the different disorder strengths with the 
parameters L = IQ, N = 60, and (f) = 0.007. 

FIG. 4: (Color online) The conductance G vs. the energy El for different width N at the moderate 
disorder strength W = 2t and L = 16. The panels (a) and (b) are for the magnetic field strength 
= and = 0.007, respectively. 

FIG. 5: (Color online) The conductance G vs. the width N of the graphene nanoribbon with 
El = O.lt (solid curve), 0.2t (dashed curve), and 0.25t (dotted curve). The panels (a) and (b) 
are for the magnetic field strengths = and (p = 0.007, respectively. The other parameters are 
W = 2andL = 16. 

FIG. 6: (Color online) The conductance G vs. the energy El for different width N. The parameters 
are L = 12, the disorder strength = [in (a) and (c)] and W = 2t [in (b) and (d)], and the 
magnetic field strength = [in (a) and (b)] and = 0.007 [in (c) and (d)]. 

FIG. 7: (Color online) The conductance G vs. the width N of the graphene nanoribbon for 
the different energy El- The panels (a) and (b) are for the magnetic field strengths = and 
= 0.007, respectively. The other parameters are W = 2t and L = 12. 

FIG. 8: (Color online) The normal tunneling coefficient T(a;) [in (b) and (d)] and Andreev refiection 
coefficient Ta{uj) [in (a) and (c)] vs. w for different on-site energy El- The other parameters are 
N = 50, W = 0, and the magnetic field strength = [in (a) and (b)] and = 0.007 [in (c) and 
(d)]. 

FIG. 9: (Color online) The normal tunneling coefficient T(lu) [in (b) and (d)] and Andreev refiection 
coefficient Ta{uj) [in (a) and (c)] vs. u for different on-site energy El- The other parameters are 
N = 50, L = 16, W = 2t, and the magnetic field strength = [in (a) and (b)] and = 0.007 [in 
(c) and (d)]. 
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